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This  paper  presents  a general  theorem  on  the  invariant  behavior  of 
a certain  function  of  a matrix.  It  then  shows  the  importance  of  this 
result  principally  by  using  it  to  derive  properties  of  certain  maximum 
likelihood  estimates  which  arise  idien  considering  problems  such  as  the 
location  of  a moving  object  being  surveyed  from  a moving  observatory 
when  all  data  on  location  are  subject  to  stochastic  error.  This  problem 
is  important  in  tracking  objects  either  from  an  observatory  satellite  or 
from  a transport  plane  beeiring  ground  seeking  radar.  Some  applications  to 
this  situation  are  made. 
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1.  INTRODUCTION 


This  paper  stems,  in  part,  from  two  earlier  proprietary  reports,  viz. 
Saunders  (I965)  and  Saunders  and  Johnson  (1964),  dealing  with  statistical 
problems  arising  from  the  estimation  of  the  position,  heading  and  velocity 
of  a moving  object  using  data  which  are  subject  to  statistical  error. 

These  data  were  presumed  to  have  been  obtained  from  observations  made 
during  one  overflight  -vdien  the  exact  location  of  the  observatory  platform 
is  not  precisely  known  with  respect  to  the  ground.  The  solution  of  this 
problem  was  originally  intended  to  assist  in  the  determination  of  both 
the  travel  and  location  of  ships  \dien  using  data  obtained  from  an  observa- 
tory satellite.  Recently  the  same  mathematical  problems  have  arisen  with 
the  introduction  of  ground  mapping  radar  vhich  is  being  born  by  airplanes 
and  used  in  the  observation  and  prediction  of  position  of  moving  land 
vdaicles. 

In  the  first  section  a theorem  on  the  b^iavior  of  a certain  function 
of  a matrix  is  stated  and  proved.  In  the  following  sections  a simpi.e 
model  with  normal  errors  for  the  moving  traget  from  a moving  observatory 
is  given  and  the  maximum  likelihood  estimates  of  target  position  are 
obtained.  The  theorem  is  then  utilized  to  yield  certain  invariance 
properties  of  these  estimates. 


2.  THE  GENERAL  THEOREM 

We  now  state  and  pi'ove  a general  res\ilt  on  the  behavior  of  a particular 
function  of  positive  definite  matrices. 
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(2.1) 


independent  cf  t.  If  the  nvuiierator  is  independent  of  t for  all  u 
then  the  denominator  will  also  be  independent  of  t by  taking  n = 0. 

Using  the  fact  tliat  if  A is  non-sing-jlar 

IJ,  l\  - IaIId-b-a'^b!  , 

we  obtain  for  the  numerator  of  (2.])} 

I Zq+u' ul  • i (t+c  ) hi’ u - 

F ec.u« u+t  u ) ] u ) “^  + (t-K:  )a' u)  | . 

The  first  term  is  independent  of  t.  The  matrix  in  the  second  term  is 

£2+2tZi+t^ZQ  + (t+cfu’u  - {i:^+cu'u}(Z^+u'uf\z^+cu’u) 

- 2t  (Z-j^+cu' u)  - t''(Z^-ni'u)  , 

viiich  upon  expandon  and  simplification  reduces  to 

2 -1 
Zg-i-c  u’u  - (Z^+cu' u)  (Z^+u’ u)  (z^+cu'u)  . 

This  expression  is  independent  cf  t,  which  completes  the  proof.  ll 

3.  Tlffi  MOTEL  AND  ITS  ASSUMPl'IONS 

We  shall  speak  of  the  object  (or  target)  and  the  observatory,  including 
thereby  all  applications,  vritli  the  understanding  that  both  points  are  moving 

3 


with  respect  to  a given  coordinate  axis.  Oiu”  first  task  is  to  derive  the 


appropriate  density  of  the  observations  of  target  positions  relative  to  a 
fixed  coordinate  axis,  as  determined  from  the  observatory. 

We  now  specify  precisely  the  assumptions  on  vdiich  our  analysis  is 
based.  These  are: 

1°  The  observatory  moves  at  constant  hei^t  linearly  above  the  plane 
with  a known  constant  velocity. 

2°  Tlie  object  moves  linearly  in  the  plane  at  a constant  but  unknown 
velocity. 

3°  The  estimated  position  coordinates  of  the  observatory  over  the 
plane,  as  determined  from  the  ground  at  a given  time,  are 
bivariate  normal,  random  variables  with  known  covariance. 

^4-°  The  estimated  position  coordinates  of  the  object  in  the  plane, 
as  determined  relative  to  the  true  position  of  the  observatory, 
sire  bivariate  normal  random  variables  and  successive  observations  of 
such  relative  positions  are  independent. 

3°  Time  between  successive  observations  can  be  measured  with  sufficient 
accuracy  so  that  errors  of  position  due  to  time  inaocuracy  are 
negligible. 

6®  The  parameters  of  the  covariance  matrix  of  the  observations  of 
I object  position  relative  to  the  true  observatory  position  can 

be  determined  from  bearing  angle  and  range  data. 

On  a single  overflight,  an  observatory  may  make  several  observations 
of  the  position  of  an  object  and  our  first  problem  is  to  determine  the 
joint  distribution  of  the  observations  of  the  target  position,  using  the 
fact  that  the  target  position  relative  to  the  observatory  is  subject  to 
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observutional  error,  as  is  the  estinate  of  the  observatory  positions 
relative  to  the  ground.  Ihus,  at  a given  tiae  t^  we  assume  the 
ohservatorj’  is  at  some  position,  say  P^,  and  the  target  at  some 
position  (.u  both  in  a plane  located  with  respect  to  a given  coor- 
dinate system.  However,  the  position  of  the  object  an  observed  by 
radar  from  the  observatorjr  is  subject  to  chance  error  and  hence,  the 
radar  estimate  of  the  object  position  from  the  observatory  position  is 
a random  variable,  say  X^.  Now  a radai’  measiirement  from  the  ground 
at  time  of  ttie  observator;/  position  on  the  given  coordinate 

system  is  also  a random  variable,  call  it  Y. 

From  assumptions  1°  and  2°  we  have  that  the  observatory  follows  a 

linear  path  in  the  plane,  say  P,  = ; + et,  as  does  the  object,  say, 

*c 

^ =:  a + Pt.  (Greek  letters  denote  points  in  the  plane.) 

Vathout  loss  of  generality,  we  can  select  our  coordinates  so  that 
the  first  coordinate  of  e is  in  the  direction  of  observatory  travel 
and  hence,  the  second  coordinate  of  e is  zero.  Moreover,  by  assumption 
1“ , the  first  coordinate  of  e is  known. 

Again,  following  the  general  mathematical  assumptions  3'*  and  4®, 
we  have  that  X^,  for  i = 1, .,.,n,  and  Y are  bivariate  normal. 

More  specifically,  X^  has  mean  vector  and  known  covariance 

matrix  C7^,  whereas  Y han  mean  vector  p„  and  covariance  matrix  D~^. 
That  is,  for  i = 1, ...,n 

~ TI(H^-P^,C"^)  and  Y ~ TiCOq^d”^)  . 

The  data  obtained  from  the  observatory  yield  the  observations 

Z.  = X. +Y  for  i = l,...,n  and  in  this  section  we  seek  to  derive  the 
11 

joint  density,  f,  of  Z=  (Z^^, . . . ,Zj^).  Since  Z^  = it  follows 
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» lU  W iJBii 


that 


and 


> 


Var(Z^)  = +d'^ 

i, j = 1,. . . ,n.  Hence 


Cov(Z^,Z  ) = d"^  for  i j 


Z = Cov(Z)  = . 


vAiere  8^^  is  tlie  Kronecker  delta.  Setting 


it  follows  that 


0.1) 


''ij  = Vi- 


This  inversion  is  a consequence  of  the 
Lemma:  If  are  non-singular  symmetric  matrices  of  the  same 

dimension  and 


H — diag  • j Qj^)  + . 


^ 1 ) j 


then 


h'^  = diag(Q'\...,0'^)  . (I  + Z OTV 

• , d.  ~1 

)«i:X 


See  Householder  (1964),  p.  124. 
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Specifically,  tlie  choice  of  = C~^,  = d"^  yields  the  result 

claimed. 

Jfoti;^  that  follows  that  A is  itself  symmetric. 

Except  for  elementary  details  we  have  ootained  the  basic 

Theorem  2:  'Hie  distribution  of  the  observations  of  target  position 
Z = (Zj^,. . . /,  obtained  from  the  observatory  is 


Z ~ n.(v,E)  , 


with  density  function 


(3.2) 


’■'(z)  = exp[-  ?r(z-v)A{;,--v)’ ] , 

(2nf- 


where 


A = Z is  defined  by  (3.1)  , and 


(3.2) 


IaI  ^ 


bi  n !c. 


D + Z C. 


'Iheorem  2 tells  us  that  a single  observation  of  observatory  position 
from  a bivariate  nonnal  distribution  relative  to  the  ground,  combined  with 
n observations  of  the  object  position,  vhich  are  normally  distributed 
relative  to  the  observatory  position  jdelds  a Joint  normal  distribution  of 
observed  target  positions  relative  to  the  ground.  Our  next  problem  is  to 
find  a best  estimate  of  the  target  coirrse  using  this  normal  distribution 
of  error  and  the  assumption  that  the  target  is  moving  linearly  at  a constant 


k.  THE  MAXIMUM  LIKELIHOOD  ESTIMATES  AND  THEIR  DISTRIBUTION 


Thus,  we  want  to  find  the  maximum  likelihood  estimate  (MLE),  ja  , 
of  future  target  positions  as  a function  of  t.  By  properties 

of  the  MLE, 

A A 

= a + &t  , 

A A 

vdiere  O'  and  P are  the  MLE's  of  o and  P.  To  complete  the  picture 

A A 

we  need  to  know  the  distribution  of 

Recall  that  for  i = l,...,n. 

If  we  define  y = a +■  et^,  k = g-E  and  obtain  the  MLE's  of  y and  'C, 
then  we  can  easily  obtain  the  MLE's  of  o:,g. 

The  likelihood  function  is,  with  ^ for  i = 1, ...,n, 

L = log)El  - n log(2jr)  - ^ f.A.  . . 

i,J^  ^ ^ 

Introduce  the  notation  for  the  coordinates  of 

then  writing  = ^S2k^^2k^  ^ where  5^^  is  the  Kronecker 

delta  we  obtain,  by  taking  partial  derivatives,  the  two  equations  for  k= 1,2, 


^ J 


0'.2) 


^ Z (t.  LA.  .8'  + t.8,  A.  ./•'.)  . 
2 ' 0 T.  ij  k 1 k 13  'j'  * 
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Realizing  the  second  term  in  each  expression  is  a scalar,  and  tiius  equal 
to  its  own  transpose,  the  two  equations  in  (4.1 ) and  the  two  equations 
in  (4.2)  can  be  written  matrix  form,  upon  equating  to  zero,  as 


(4.3) 


(4.4) 


Note  that 


3,  ' 

i ’ J 


^.+t.  = z.  + z . - 27  - K(t.  +t.)  , 
■1  -J  1 J ' 1 ’ 


t.t.+t.t.  =t.z.  +t.Js.  -7(t.+t.)  -2/ft.t.  . 
0 -a  1 3 X 1 ,3  13  13 


By  substitution,  we  obtain  from  (4.3)  and  (4.4)  the  equations 


(4.5)  S (z  +z  )a..  = 27  S a.  + /C  Z (t.+t  )a.  , 

i,3  ^ J iJ  i,3  iJ  i,j  1 d 13 


(4.6)  Z (t  z +t  z )a.  = 7 Z (t  +t  )a.  . + 2;f  z fept  a.  . 

i,j  3 1 1 3 13  i,j  ^ J i,3  1 3 13 


Writing  these  in  matrix  notation,  we  have 


(4.7) 


^31  ^12  \ 

(r,x)  . 

» ^21  22 1 -L 
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or  in  matrix  notation,  with  the  obvious  definition, 

^ = ZU  . 


By  Theorem  2,  Z is  U(v,  a"^)  and  t is  U(vU,  U'a"^U).  Since 
(7,<)  = >|fS,  (7,/c)  is  n(vUS,Q),  i*,ere  Q = S'U'a'^US.  We  now  prove  that 

VUS  = (7,«r)  . 

To  do  this  we  maice  use  ofttie  identities  which  follow  from  Z~^S  = I. 

The  first  component  of  the  vector  vUS  is 


1^1 


^^^11^11  ^12^21^  '^^^21^11  ^22®21^  ~ ^ 


The  second  component  of  the  ’.'•ector  vUS  is 


(r^t, 

J i^ J 


""  ’'^^11^12  ^12^22^  '^^^21^12  ^22®22^  - * 


which  completes  the  proof  of  (4-.  11), 


11 


J 


We  have  Just  obtained 


Theorem _3_.  The  MLE,  of  (0£,P)  has  components  given  by 

^ A A 

« = 7-etQ  , p = ^f+e  , 

where  {y ,<)  are  defined  by  (1.. 9).  Furthermore  (5,P)  ~n((a,P),Q)  vith 

Q = S'U'a'^US  = I ^ 

\^1  ^2/ 


From  tliis  we  obtain 

Corollaxy  1;  if,  from  conditions  of  symmetiy,  the  additional  assumption 

that 


A'a"^'  =a 

holds,  then 

Q = S'lnA’^US  = s . 

Proof.  From  (4.10) 


U = 


'■'ll  - ''in 


•••  F 

nl  rm 


I t^I 


Bince  F.  . F,  it  follows  that  Fa'^F  = F,  and  hence 

2 . 2 t F \ 


U'A  "U  = J 


/F.  . 

t F 

J"iJ 

Waj 

t.t.F. . 
1 J ij 

= S 


2 t.F  z t,t.F_ 

^ i,J  ^ J 


' = S.  I 
12 


The  assertion  follows  since  S 


Corollary  2:  The  MLE  of  the  true  target  position  at  time 


t is 


vfhere 


A A A 

= a + pt 


Applying  Theorem  1,  we  have  obtained  a result  of  great  practical 
importance,  namely. 

Theorem  4;  The  estimate  p^.  and  its  covariance  matrix  B.  are 

— — — - X#  ti 

invariant  under  location  and  scale  change  in  time. 


5.  SOME  SPECIALIZATIONS  OP  PRACTICAL  IMPORTANCE 


Assume  that  a given  time  t^  the  observatory  is  at  point  and 
traveling  in  a strai^t  line  at  a known  constant  velocity.  It  observes, 
at  a true  bearing  angle  6^  and  range  R^,  the  object  which  is  located 
at  position  p^  obtaining  the  random  variable  ~ ]X(p^-p^,c”^),  See 
Figure  1. 


direction  of 
travel  at  known 
velocity 
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From  elementary  geometrical  argiunents,  as  in  Saunders  and  Johnson  (196*4- ), 
which  need  not  lie  given  here,  we  have 


(^.1) 


P P 2 2 

! CT  cos  0.  +CT..sin  0. 
/ Ri  1 Ai  1 


\ '^Ri"'^Ai  . „ , 

- r; sin2  0. 

(.  1 


2 2 
On.  -a,. 

^ sin  2 0^ 


2 2 2 2 
Oj^sin  0^  +cT^cos  0^ 


cr 


v±Lere  at  time  fc^,  cr^  is  the  standard  deviation  of  the  range  error, 
is  the  standard  deviation  of  the  azimuth  error  and  both  are  known 
functions  of  the  range  Ri,  all  in  accord  with  Assumption  6®. 

We  also  assume  that  the  matrix  D is  diagonal  and  known. 

If  it  is  true  that  the  distribution  of  radar  errors  is  constant  in 
time,  that  neither  the  direction  of  travel  of  the  observatory  platform 
nor  its  position  relative  to  the  target  will  influence  the  covariance 
matrix  of  observations,  then  C.  = C for  i = l,...,n  and  A. . = A... 

^ a j ji 

Qy  a strcu^tforward  calculation,  with  the  time  chosen  so  that  t"  = Et^/n  = 0, 
we  obtain 


= nC-n^C(D-tnC)"^C  , E = ^21  = ° ’ 202  = n t^  C 


'*^1  " ^^11  ’ ^2  = n tz  C , 


(where 


- E z^,  “ n ^ ^i^i’  ^ “ 2:t^/n  ) and  hence 


^11  - n ° , 3,5  - So,  - 0 , S, 


12  “ “21 


^22 


,-l 


n t 


1*+ 


FiDm  (^.9)>  the  maximum  likelihood  estimates  are  given  by 


A _ A 


7 = z,/f=tz/t  , 


and  hence  by  Uieorem  3 and  Corollary  it  follows  that  \i  - a r $t,  where 

Xi 


A — 


a = z - et^  , ^ ^ + e . 


Hie  estimator  is  unbiased  and  has  covariance  matrix 


t^ 


In  the  case  that  C is  diagonal,  the  matrix  B is  also  diagonal. 

t 

In  the  circumstance  that  the  observations  are  symmetrically  spaced 
in  time  and  the  origin  is  chosen  so  that 


ti  + *n+i-i  0 ^ + ®n+l-i  “ " ' i - l,...,n  , 


it  follows  by  (5.1)  that  ^ is  a diagonal  matrix;  the  matrices 
n 

D + Z C^,  are  also  diagonal,  liiereas  Z^j^  is  contradi agonal. 
This  facilitates  the  computation  of  S'U'AUS  which  will  not  in  general 
reduce  to  S as  in  Corollary  1. 


6.  ESTIMATION  AND  CONFIDENCE  INTERVALS  FOR  SPEED  AND  HEADING 
The  parameters  of  interest  when  tracking  a moving  object  etre  the  speed 
and  heading,  and  in  this  section  we  consider  the  accuracy  with  vAiich  they 
can  be  estimated.  The  estimate  of  the  true  position  at  any  time  t. 


" r .. 


A A A 


(t)  =a  + t<t  = (inj^(t),  > 


\d)ich  we  write  in  coordinate  form,  assuming  it  follows  a linear  path. 
The  true  velocity  s,  heading  angle  h,  are  given  by 


s = {[m^(t)l  + [m^(t)3  ) , h = arctan 


These  relationships  are  seen  to  be  simply;  if  ^ = (b^,b^),  then 


(7.1) 


b^  = s cos  h , b^  = s sin  h . 


Thus  by  analogy  we  have  the  equations 


(7.2) 


b^  = V cos  <5  , bg  = V sin  $ 


defining  the  random  variables  V (for  velocity)  and  $ for  the  heading 
angle,  which  are  estimates  of  the  true  speed  s and  the  true  heading  h. 

Theorem  5,  it  follows  by  known  results  on  the  marginal  distribution 

A 

of  normal  variates  that  & ~ joint  density  of  V is 

found  by  simi^ly  making  a transformation  to  polar  coordinates.  Ihis  yields 


(7.=^) 


0 < V < 00 


g(T,vlh,s)  = / exp(-  I r,  0:2  (;•)  > 

2A\\  ^ ^ ‘ 0 < cp  < 2;; 


where  ^ = (v  cos  rp  - s cos  h,  v sin  q)  - s sin  h).  'This  density  can  be 
used  to  study  the  distribution  of  velocity  and  heading  estimates  that 
could  arise  under  infrequent  headings  and/or  hi^  velocity. 


! 

I 


t 


I 


one  mitfjit  desire  senarate  conr'idence  intervals  on  the  heeuding 
and  on  the  velocity.  However,  if  we  proceed  to  find  the  marginal 
densities  of  v and  '(i  from  (7.'')  we  see  that  each  density  has  both 
parameters  h and  s, 

Tims  a confidence  interval  for  the  velocity  can  be  constructed 
only  if  we  know  the  trae  heading  h.  Likewise  a confidence  interval 
can  be  found  for  the  tnae  headir^  if  we  know  the  true  velocity. 
presence  of  the  nuisance  parameters  prevents  us  from  obtaining  confi- 
dence intervals  separately  when  both  parameters  are  unknown. 

However,  we  can  obtain  a joint  confidence  region  for  ^h,s),  vhich 
is  somewhat  inefficient,  as  follows j irom  well-known  results  on 
the  chi-square  distribution  of  the  quadratic  form  of  normal  variates 
we  have 

(7.^0  < x;(p)}  1-p  , : 

c cl 

2 

where  X^Cp)  can  be  easily  calculated  for  any  0 < p < 1, 

The  (random)  elliptical  region  W,  so  defined,  determines  a 

100 (l-p)  percent  confidence  region  for  P.  We  seek  the  smallest  area 

in  polar  coordinates  which  is  the  Cartesian  product  of  intervals  and 

■X- 

contains  Wj  call  it  W . See  Figure  2. 

1 

A Figure  2 I 

^ I 


Joint  confidence  intervals  for  heading  angle  C 
and  velocity  V, 


La 
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Witii  the  obvious  notation  for  the  maximum  and  minimum  of  argument 
and  modulus)  we  write 


(7.!?) 


W ( (0,r):  < 0 < <5^,  < r < V^}  . 


rfow  (r  cos  t',  r sin  tO  e W implies  (e^r)  e W*,  P((h,s)  e W*)  > 1-p  . 

'Hie  task  to  which  we  now  address  ourselves  is  the  determination  of 
a functional  representation  for  the  random  variables  (iJ).,  V.  ),  i = 1,2. 

If  we  denote  the  elements  of  the  symmetric  matrix 


’^o(p) 


*^11  ^^12  ! 


h2  ^^2  / 


i ^ 


then  we  can  write  from  (7.^4-) 


(7.6)  'lll-^"'^l)^  + 2 q^2(x-b^)(y-b2)  + ^ 


We  want  to  find  the  maximum  and  minimum  of  the  functions 


(7.7) 


f(x,y)  = y/x  , g^x,y)=x^+y^ 


subject  to  the  restriction  (7.6). 

Proceeding  by  analytic  geometry  or  by  the  method  of  Lagrange  multi- 
pliers, leads  to  the  solution  of  a quartic  equation.  A simpler  method 
seera-s  to  be  the  following:  transform  (7.6)  to  the  variables  (u,v)  via 


X = bj^  + u cos  0Q  - V sin  G , 


y = b^  + u sin  Oq  + v cos  G^  , 


where 
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-q 


tan  20„  = 


12 


0 q22“'iii 


Thus  we  want  to  extremize  the  functions  of  (u,v)  obtained  by  substituting 

2 2 

into  (7*7)  subject  to  the  restriction  &^u  + b^v  = 1.  Elementary  analysis 
shows  that 


2 2 
^1  " '^11  ^0  *^12  ^2  ^ 


&2  = ®o  " *^12  *^2  ®0  * 


Now  let 


and  set 


cos 


qp  = u , sin  cp  = /6^  V , 


cos  6 


0 


sin  B 


0 


, i - 1,2  . 


The  problem  then  becomes  that  of  extremizing  the  functions 


It  cos  cp  + t>T  sin  qp  + b„ 
f(<P)  - ^ 7^ 


COS  qp  - sin  cp  + b^ 


and 


A 2 A 2 

e(l>)  = (r^cos  cp  - Ig  9 + COS  <p  + sin  (p  + b2) 

for  the  range  of  values  0 < (p  < 2x,  It  is  straightforward  to  show  that 
f (qp)  =0  if  and  only  if 
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(7.B)  i sin  cp  + cos  cp  = 0 , 

viiere 

'^11  " '2^1  ^ ’b.^  " "^1  " ^2^1  ■ ^1^1  » '^2  ""  ^2^2 

with  capital  letters  denoting  those  random  variables  which  are  functions 
of  the  estimates. 

Now  let  be  the  unique  random  angle  for  vAiich  sin  = T^/t, 

cos  Oq  - T^/T,  where  T = . Then  (7«8)  becomes  simply,  by 

using  the  trigonometric  identity  for  cosine  of  a sum,  cos  = -6j^/t. 

Letting  be  the  angle  on  (0,3t)  such  that  0 = arc  cos(-&j^/t) 
we  have,  since  cosine  is  an  even  function,  ^2  ~ ^ 

trie  two  solutions  of  '7.3).  Thus  we  have  as  the  limits  on  the  true 
bearing  angle  as  in  (7.5),  the  angles 

<5^  ^ arctan  f(S2j^)  , = arctan  f(a^)  . 

This  acccmplishes  the  location  of  local  extremum  of  f.  Graphical  examination 
of  the  function  or.  C0,2n)  may  be  necessary  to  see  if  these  are  the  true 
extrema. 

The  exurcmum  of  g is  obtained  in  a manner  similar  to  that  for  f. 
ihis  yields  g' (v)  =0  if  and  only  if 

'7.9)  sin  2T  + cos  2cp  + 1^  sin  (p  + cos  cp  = 0 

'/there 
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^2  ^1^2  ~ "'^2 


2 2 2 2 
^6,  = - ^2’ 


A 

L^  =b 


1 '1  ' ^2^1  ’ 


1^2  ^1^'2  ~ ^2^2  ’ 


and  again  upper  case  letters  denote  random  variables.  The  solutions  of 
(7-9),  call  them  in  general  be  found  explicitly,  but 

must  be  found  numerically.  This  is  easily  carried  out. 

Having  found  the  solutions  , i - 1,2  the  limits  on  the  velocity, 
vdiich  were  sou^t,  become 


= miniVg(r^  ) , (/iTTp")  - max (i/g (f^  ) , /g(r^))  . 

A circumstance  of  practical  interest  occurs  \dienever  the  observations 
are  drawn  symetrically  on  each  overflif^t.  In  this  case  0^  =0.  ^■|_  = 9^1' 

^2  • <>22'  ^12  “ 


^ ^2  ~ 9 — 


i/9 


11 


^^1  " " *122  ^ ° ' "^2  • 


Thus  the  equation  (7.9)  reduces  to 


(7.10) 


>v  A 

11  1 ^1  ^2 
- (;^ — ) Sin  2cp  sin  cp  = — — cos  (p  . 

^2 


^ *^11  '^22 


If  we  further  specialize  and  assume  that  q^^^^  = q^p  we  see  from 


(7.10)  that  the  equation  to  be  solved  is  merely  tan  cp  = bg/b^ 
which  has  solutions 


A 

COS  fp  = i > sin  cp  = + , 
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A 

^2 

- V 


From  the  defini- 


where  V is  the  estimate  of  velocity  defined  in  (7.2). 
tion  of  g previously  given  we  see  that  implies  that  the 

function  we  were  extremizing  was  merely 


cos  qj 


\2 


sin  cp  '' 




and  we  obtain,  directly  from  the  definition  of  and  previously 

made,  the  simple  answer 


V = V + = V ^ 

v/q-,1  v/q^^ 
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